Fluctuational internal Josephson effect in topological insulator film 
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Tunneling between opposite surfaces of topological insulator thin film populated by electrons and 
holes is considered. We predict considerable enhancement of tunneling conductivity by Cooper 
electron-hole pair fluctuations that are precursor of their Cooper pairing. Cooper pair fluctuations 
lead to critical behavior of tunneling conductivity in vicinity of critical temperature with critical 
index v = 2. If the pairing is suppressed by disorder the behavior of tunneling conductivity in vicinity 
of quantum phase transition is also critical with the index n = 2. The effect can be interpreted as 
fluctuational internal Josephson effect. 
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I. INTRODUCTION 

Cooper pairing of spatially separated electrons and 
holes was predicted in the system of semiconductor quan- 
tum wells more then thirty years agoi. Later it was 
observed in quantum Hall bilayer at total filling factor 
vt = I that can be presented as the system of spatially 
separated composite electrons and composite holes {seem 
and references therein). After graphenc discovery Cooper 
pairing of Dirac electrons and holes in the structure of in- 
dependently gated graphenc layers has been proposed^"—. 
Recently possibility of Cooper pairing of Dirac electrons 
and holes was predicted^ii in thin film of topological in- 
sulator, new unique class of solids that has topologically 
protected Dirac surface states^. The electron-hole pair- 
ing in that system is the realization of topological su- 
perfluidity and hosts Majorana fermions on edges and 
vorticesiiiS that is the topic of extraordinary interes t^^'^^ . 
Also Cooper pairing can lead to a number of interesting 
physical effects including superfluidity^, anomalous drag 
effeclii^, nonlocal Andreev reflectioni^ and prominent in- 
ternal Josephson effecti^. 

Internal Josephson effect manifests itself in colossal en- 
hancement of tunneling conductivity at zero voltage bias 
and it has been observed in quantum Hall bilaye r^^i^^ . It 
is caused by coherence between electron and hole states 
which leads to a coherent tunnel current depending on 
the phase of electron-hole condensate. Microscopical 
and macroscopical description of internal Josephson ef- 
fect was addressed in a number of interesting theoretical 
papersi^"— . 

Cooper electron-hole pairs can appear above critical 
temperature as thermodynamic fluctuations. They lead 
to logarithmic divergence of a drag conductivity as func- 
tion of a temperature^^i^ and a pscudogap formation in 
single-particle density of states of electrons and holes'^-. 
Manifestations of Cooper pair fluctuations in tunneling 
have not been consider previously. The pairing fluctua- 
tions can transfer coherent tunnel current and enhance 
tunneling conductivity above critical temperature. Since 
amplitude of pairing fluctuations increases in vicinity of 
critical temperature, as fluctuations of ordered state do 



for different phase transitions, one can expect the sig- 
nificant enhancement of tunneling conductivity. The de- 
scribed effect can be called fluctuational internal Joseph- 
son effect and in this work we develop its microscopic 
tlieory. This effect is considered in topological insula- 
tor thin film and its peculiarities in other realizations of 
electron-hole bilayer are discussed in Conclusions. 

The rest of the paper is organized as follows. In Section 
2 wc briefiy discuss the model and microscopical descrip- 
tion of Cooper pair fluctuations. In Section 3 tunnel- 
ing conductivity of topological insulator fllm is calculated 
and Section 4 is devoted to conclusions. 



II. COOPER PAIR FLUCTUATIONS 

We consider topological insulator thin film with in- 
dependently gated surfaces occupied by excess electrons 
and holes. In the single-band approximation the Hamil- 
tonian of the system is given by^ 
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p is annihilation operator for a electron on the 



Here a 

surface with excess of electrons and b 



is annihilation 



operator for a electron on the surface with excess of 
holes^^; ^p = vpp — E-p is Dirac dispersion law in 
which v-p and Ep are velocity and Fermi energy of 
electrons and holes. Wc consider the balanced case 
since it is favorable for Cooper pairing and the pair- 
ing is sensitive to concentration mismatch of electrons 
and holes. U{q) is screened Coulomb interaction be- 
tween electrons and holes (see^ for its explicit value) and 
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/+q/2) is angle factor 



that comes from the overlap of spinor wave functions of 
two-dimensional Dirac fermions. 

For microscopical description of Cooper pair fluctua- 
tions we introduce Cooper propagator. It corresponds to 
the two-particle vertex function in the Cooper channel^! 



and satisfies the Bethe-Salpeter equation depicted on Fig. 
1 (a). In Bardeen-Cooper-Schrieffer (ECS) approxima- 
tion its solution can be presented in the form 
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where U' is the Coulomb coupling constant^ and nj?(aj) 
corresponds to electron-hole bubble diagram that can be 
interpreted as Cooper susceptibility of the system. Direct 
calculation of this diagram leads to 
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Here 7 = {ja + lb)/'2 is disorder caused Copper pair 
decay rate equal to half-sum of phenomenological intro- 
duced decay rates of electrons and holes; Tq is the critical 
temperature of electron-hole pairing without disorder; fp 
is the density of states of electrons and holes on the Fermi 
level; ^{x) is the digamma function. In the absence of 
disorder Tdio) = at the critical temperature Tq indi- 
cating Cooper instability of the system against Cooper 
pairing. Critical temperature for disordered system T^ 
satisfies the following equation 
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(4) 



This equation has nontrivial solution if 7 < 7c, where 
7c = 0.89To is the critical Cooper pair decay value. In 
opposite case the pairing is suppressed by disorder. The 
value 7c corresponds to quantum critical point at zero 
temperature. 

In a vicinity of the critical temperature T^ the expres- 
sion for Cooper pair propagator ([3]) can be approximated 
in following way 



r^(^) 



1 



4^Td 



t-F 47r(T - Td) - iuj^'{4 



2-kTa 



) 



(5) 



In a vicinity of quantum phase transition at zero temper- 
ature Cooper pair propagator is given by 
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III. TUNNELING CONDUCTIVITY 

The tunneling of electrons between the opposite sur- 
faces of topological insulator thin film with conserving 
momentum can be described by the following Hamilto- 
nian 



i/T = T + r+ = ^ (f6+ap + ra+&p) , 



(7) 



where t is the tunneling amplitude. For a calculation of 
the tunneling conductivity we use linear response the- 
ory in which the tunneling conductivity i7t(1^) at finite 





T 



FIG. 1. a) Diagrammatic representation of the Betlie-Salpeter 
equation for the Cooper propogator Fc; b) Feynman diagrams 
for a tunnel conductivity. Solid (daslied) line corresponds to 
electrons on the surface of TI film with excess of electrons 
(holes). 



voltage bias V can be presented in the form of Kubo 
formula^S 
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where the retarded response function Xt('^) ^^'^ be ob- 
tained by analytical continuation iri„ — > w + i(5 of 
X^(il7„) that is given by 



X^i^^n) = -^J_ dr 



•^(rMT(T)T+(0)). (9) 



Here Tm is time-ordering symbol for imaginary time r 
and rin = 2TrnT is bosonic Matsubara frequency. In the 
system of noninteracting electrons and holes the x^(ir2„) 
corresponds to the first diagram on the Fig 1 (b) leading 
to Xt('^) = I^Pn^C'*^)- Hence the tunneling conductivity 
for noninteracting electrons and holes ctto is given by 
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Its value CT^Q(7,r) at zero bias is given by 
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and at low temperatures T <C 7 it transforms to 
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Introduction of the electron-hole Coulomb interaction in 
the ladder approximation leads to the additional term 
for tunneling conductivity that corresponds to the sec- 
ond diagram on the Fig 2 (b). Infinite series of electron- 
hole scattering diagrams in Cooper channel correspond 
to Cooper propagator hence the second term can be at- 
tributed to the pairing fiuctuations. The tunneling con- 
ductivity for interacting electrons and holes cttc is given 

by 




V (mV) 



FIG. 2. (Color online) Tunneling conductivity ctt as a func- 
tion of bias voltage V for noninteracting (dashed line) and 
interacting (solid line) electrons and holes for 7 = 0.2 K and 
r = 0.2K. 
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The denominator of P^ coincides with that in the 
Cooper propagator ^ and tends to zero in a vicinity of 
critical temperature T^- Hence the Cooper pair fluctua- 
tions lead to critical behavior of tunneling conductivity 
in vicinity of critical temperature and quantum critical 
point. In vicinity of critical temperature Tj tunneling 
conductivity at zero bias is given by 



f^Tc =crTo(7:^d 
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In vicinity of critical temperature it diverges as (Tt ~ 
(T — Td)^'^ with the critical index :/ = 2. At zero tem- 
perature in the vicinity of the quantum phase transition 
at 7 = 7c the tunneling conductivity is given by 
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It diverges as (Tt '^ (7 — 7c)~'' with critical index fj, = 2. 
The formulas p4|) and p^ arc the main result of the 
article 

For numerical calculation of the tunneling conductiv- 
ity at finite voltage bias we use formulas ([5)1. P^ and the 
following set of the parameters^ Tq = O.IK, /j, = 5 meV, 
t = 2 fieV. The phase diagram of the system is pre- 
sented on the inset of Fig. 3. If Cooper pair decay rate 
exceeds value 7c = 0.09 K then the electron-hole pairing 
is suppressed by disorder. 

Calculated tunneling conductivity for 7 ~ 0.2 K and 
T = 0.2 K is presented on Fig. 2. Due to restrictions con- 
nected with energy and momentum conservation the de- 
pendence has prominent peak. Such peak was predicted 
and observed also in electron-electron bilaye r^^i'^^ and it 
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FIG. 3. (Color online) Half-width of the tunnel conductivity 
peak Vhw for interacting electrons and holes as function of 
temperature T for different values of Cooper pair decay rate 
7. Inset: The phase diagram of the system in which BCS 
denotes the pairing state. The arrows on the phase diagram 
correspond to the dependencies on Fig. 3 and Fig. 4. 



is peculiarity of tunneling between two two-dimensional 
systems. Coulomb interaction between electrons and 
holes considerably enhances tunneling conductivity but 
does not change qualitatively its dependence on external 
bias. Half- width and height of the peak for interacting 
electrons and holes are presented on Fig. 3 and Fig. 4, 
respectively. For 7 < 7c Coulomb interaction leads to 
the critical behavior of the tunneling conductivity that 
we interpret as manifestations of Cooper pair fluctua- 
tions. In vicinity of the critical temperature height of the 
peak diverges with the critical index v ^ 2 that agrees 
with the analytic results and width of the peak linearly 
tends to zero. If Cooper pair decay exceeds critical value 
7 > 7c Coulomb interaction considerably enhances the 
height and leads to reduction of the width but does not 
lead to any singularities. The peak becomes more promi- 
nent with decreasing of decay and temperature as it does 
in the model of noninteracting electrons and holes. 



IV. CONCLUSIONS 

We have shown that electron-hole Coulomb interac- 
tion considerably enhances tunneling conductivity in 
electron-hole bilayer even when the Cooper pairing is 
suppressed by disorder. The opposite situation takes 
place in electron-electron bilayer that also can be realized 
in semiconductor quantum well structure, in graphene 
double layer system and in a film of topological insulator. 
Coulomb interaction gives contribution to decay of elec- 
trons that was analyzed in^ and to additional series of 
diagrams for the tunneling conductivity. We treated the 
additional diagrams in ladder approximation (See Fig. 2- 
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FIG. 4. (Color online) The height of tunnel conductivity peak 
cr™;. for interacting electrons and holes on temperature T for 
different values of Cooper pair decay 7. 



b). If they arc omitted the tunneling conductivity at zero 
temperature and at finite bias V is given by^i^ 
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where g is degeneracy factor of electrons and tx is 
additional factor that depends on internal nature of 
electrons^. The dependence of tunnel conductivity on 
external bias contains prominent peak which becomes 
more prominent with decreasing of decay rate 7. If the 
Coulomb interaction is treated in ladder approximation 
the tunneling conductivity is given by 
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For electron-electron bilayer Coulomb interaction does 
not influences the height of the peak and leads to de- 
creasing of peaks width which is insignificant even in 
case of strong interaction vpU' ~ 1. The roles of the 
interlayer Coulomb interaction in electron-electron bi- 
layer and electron-hole bilayer are drastically different be- 
cause the correction to the tunneling conductivity of the 
primer is caused by the scattering diagrams in particle- 
antiparticle channel and the correction to the one of 
the latter is caused by the diagrams in particle-particle 
Cooper channel containing instability. 

The mean field theory we use here for description 
of fluctuational internal Josephson effect does not ac- 
count large scale fluctuation of phase of Cooper pair 
condensate. In two-dimensional superfluids phase fluc- 
tuations destroy long-range coherence and the transi- 
tion to paired state at critical temperature Tj calculated 
within mean field theory is smoothed. Moreover the 
transition to superfluid state is Berezinskii-Kosterlitz- 
Thouless transitional^ that corresponds to dissociation 



of vortex-antivortex pairs and which temperature is be- 
low Tj. Hence large scale phase fluctuations of Cooper 
pair condensate can smooth critical behavior of tunnel 
conductivity we predict here. But if the size of the sys- 
tem L is comparable with coherence length of Cooper 
pair fluctuations 4 ~ ^ivf/Tq the phase fluctuations are 
unimportant and mean field theory is well applicable. For 
Tq = 0.1 K the coherence length of Cooper pair fluctua- 
tions is of order Ic ^ 10 mkm and we conclude that the 
developed microscopical theory is applicable for samples 
of the corresponding size. 

The model we use here implies conservation of the mo- 
mentum of tunneling electron. If the momentum is not 
conserved the tunneling process creates electron-hole pair 
with nonzero total momentum of order l^ . Here It is 
character length at which tunneling matrix matrix ele- 
ment t can be considered as constant. Cooper pair is 
formed by electron and hole with opposite momenta and 
the Cooper instability is smoothed if ^c ^ ^t- For tun- 
neling between opposite surfaces of topological insulator 
thin fllm of high crystalline quality momentum conserva- 
tion can be achieved to remarkable degree. 

We have investigated the manifestations of Cooper 
electron- hole pairing fluctuations in thin fllm of topolog- 
ical insulator on tunneling between its opposite surfaces. 
The effect takes place in other realizations of electron- 
hole bilayer but each system has its own peculiarities. 

Dirac points in graphene are situated in corners of first 
Brillouin zone. Electron-hole pairing was predicted^"— in 
the system of two independently gated graphene layers 
separated by dielectric film. In that case orientations of 
the graphene lattices are uncorrelated. The distance be- 
tween Dirac points of different layers in momentum space 
is of order Uq , where oq is lattice constant of graphene. 
The tunneling of electrons between Dirac points is pos- 
sible if ^T ^ ao that corresponds to tunneling through 
impurity states or other defects. The condition Ic 3> It 
is well satisfied and the critical behavior of tunneling con- 
ductivity in double layer graphene system is considerably 
smoothed. 

Recently anomalies in drag effect in semiconductor 
double well structure that contains spatially separated 
electrons and holes were observed^i^. The analysis of 
results shows that the observed anomalies can be caused 
by electron-hole pairing not in BCS regime but rather 
in regime of BCS-BEC crossove r^^i'^^ . In that regime 
electron-hole pairing fluctuations also should increase 
tunneling conductivity of the system but quantitative 
theory of the effect is interesting and challenging prob- 
lem. The developed here microscopical theory of fluc- 
tuation internal Josephson effect is applicable for that 
system if the pairing is realized in BCS regime. 

Internal Josephson effect has been observed in quan- 
tum Hall bilayer at total occupation factor j^t = 1- But 
above critical temperature the dependence of tunneling 
conductivity on external bias does not contain any peak 
due to non Fermi-liquid behavior of composite electrons 
and holes^"— . Thus fluctuational effects in the system 



are more complicated ones and need separate investiga- 
tion. It should be noted that critical behavior of tun- 
neling conductivity has not been observed yet in experi- 
ments in that system. 

We have considered influence of electron-hole Cooper 
pair fluctuations that are precursor of Cooper pairing in 
topological insulator film on tunneling between its op- 
posite surface. Cooper pair fluctuations lead to critical 
behavior of tunneling conductivity in vicinity of critical 
temperature with critical index v = 2. If pairing is sup- 
pressed by disorder the behavior of tunneling conductiv- 
ity in vicinity of quantum critical point at zero tempera- 



ture is also critical with critical index fi = 2. The effect 
can be interpreted as fluctuational Josephson effect. 
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